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I. INTRODUCTION 

The classical Hartree-Fock (HF) approximation is a 
prototype of the modern approach to the quantum many- 
body problem related to the energy density functional 
Qtllb- When applied to complex nuclei, the density func- 
tional theory may provide a universal description across 
the nuclear chart. The pairing interaction that is present 
in nuclei as well as in fermionic condensed matter sys- 
tems is usually included in the Hartree-Fock-Bogoliubov 
(HFB) form Q. The well known deficiencies of the HFB 
approach for mesoscopic systems follow mainly from its 
non-conservation of particle number. As a result, un- 
physical features are introduced into dynamics, the su- 
perfluid phase transition appears too sharp, and the cor- 
relational energy produced by pairing might be severely 
underestimated. As was shown earlier [HIH, the pairing 
part of the problem can be solved numerically quite easily 
with the help of the seniority representation in a spheri- 
cal basis, and its exact solution significantly improves the 
results. 

It was also sketched in Q how other parts of the inter- 
action can be incorporated into the exact pairing method 
in the approximate way that reminds the HF approach. 
This can be done in an iterative fashion: the exact pair- 
ing solution using the starting single-particle basis de- 
termines the actual occupation numbers; these (in gen- 
eral, fractional) occupancies self-consistently determine, 
in the HF spirit, an improved single-particle basis where 
we again solve the pairing problem etc. until conver- 
gence. In this way both mean-field features, deformation 
and pairing, are accounted for. The main purpose of the 
current work is to further develop this Hartree-Fock plus 
pairing correlation (HFP) method that essentially is an 
intermediate step from the HF approach towards the full 
shell-model (SM) description. On one hand, we would 
want to keep the simplicity and modest computer de- 
mands as inherent properties of the mean field approach. 
On the other hand, we take into consideration pairing and 



other physical effects beyond the simple HF, or mean field 
in general. We check our approach for the scZ-shell nuclei, 
where the SM with large-scale diagonalization works per- 
fectly [f| and can serve as a searchlight illuminating the 
correct direction of motion. The success of this attempt 
will allow the extension of the approach to heavy nuclei, 
where the catastrophic growth of dimensions makes the 
complete shell model solution unrealistic. 

II. OUTLINE OF THE METHOD 

As in most mean-field approaches, we formulate this 
method as variational one. As in the shell model, we as- 
sume a general form of the two-body Hamiltonian that 
includes the single-particle term t and the (antisym- 
metrized) two-body interaction V: 

H = t ik a\a k + ^ X! Vi i kl a l a j a kai- (1) 

ik ijkl 

The variational wave function |\&) will be defined below. 
The wave function and all properties of the system follow 
from the minimization of the expectation value 

(*|H|*>. (2) 

For our test of the methods, we will take for V the USDB 
interaction from the sd-shell model @. It will allow us 
to compare the results obtained using our approximate 
method with the exact shell model calculations in the 
same single-particle model space. 

The ground state wave function |<3>) for a fixed particle 
number N can be presented as a superposition of basis 
states, 

l*) = £cyd), (3) 

where each basis state \d) is a Slater determinant which 
for N fermions can be written as usual: 

\d)=atai 2 ...at 10). (4) 
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The single-particle states </>„ can be found with the help 
of the variational principle as it is usually done in the 
HF method. The approach is actually defined by the se- 
lection of the space D spanned by the determinants \d). 
If we choose only one Slater determinant as our varia- 
tional wave function (J3]), we come to the standard HF 
approximation. If the manifold D includes all possible 
configurations, then we get the exact shell-model solu- 
tion. 

In this article our choice is determined by the pairing 
phenomenon which smears the Fermi surface and con- 
verts the Fermi-gas ground state into a superposition of 
Slater determinants. In the case of a spherical system 
with the pairing forces taken as the J = 0, T = 1 part 
of the two-body interaction {T]), we have seniority s as a 
good quantum number. For an even number of particles, 
the ground state has s = 0, while for an odd number 
s = 1. In this simple case we can construct the basis of 
Slater determinants \d) occupying single-particle levels 
\jm) by pairs, 

n *wj°>> ( 5 ) 

j;m>0 

where = (— l) J_m a]_ m is the creation operator for 
the time-conjugate single-particle state with respect to 
a ]m- Here we omit all quantum numbers except total 
angular momentum j and its projection m. 

The presence of other types of the interaction in gen- 
eral breaks spherical symmetry and brings in the de- 
formed mean field. In the case of a deformed nucleus, 
even if we had had only J = part of the two-body in- 
teraction fT} in the spherical representation, we have to 
take into consideration a broader class of pairs arising as 
a result of splitting and mixing of the original spherical 
states by deformation. Here we limit ourselves by the 
case of axially symmetric deformation, when the single- 
particle orbitals | vm) are still characterized by the angu- 
lar momentum projection m along with other quantum 
numbers v. 

According to the Kramers theorem, the orbitals \vm) 
and \v — m) are degenerate. However, the pairs may also 
be formed by the states m and — m belonging to different 
sets of remaining quantum numbers. Thus, for our basis 
Slater determinants |d) we assume the following form: 

II «L4-Jo). (6) 

i/,^;m>0 

We construct the variational wave function ([3]) as a super- 
position of the Slater determinants ^ for a given particle 
number. Using such a form we hope to correctly account 
for pairing correlations in the deformed case at the same 
time crucially reducing the dimension of space D in com- 
parison with the full shell model. Actually prescriptions 
(|5I6[) are valid only for an even number of particles. For 
an odd particle number, we use the same eq. © but 
add one additional creation operator that corresponds to 
the odd particle. The odd particle can be placed in any 



empty single-particle state, and the states are divided 
into classes with a definite value of the angular momen- 
tum projection. 

In the current application of our method we make a 
simplifying approximation treating protons and neutrons 
separately. It means that, though we consider the full 
two-body interaction ((T|) including T z — part, the vari- 
ational function ([3]) is constructed as the product of pro- 
ton and neutron parts. Clearly we are losing here proton- 
neutron correlations although their mutual contributions 
to the mean field are fully accounted for. 

The variation over amplitudes Cd with the additional 
normalization condition of the wave function, (^l^) = 1, 
leads us to the usual set of equations, 

Y,{d\H\d')C d , = EC d . (7) 

d' 

The matrix elements (d\H\d') are calculated for the de- 
terminants built on a given single-particle basis, and 
equations ((7|) are solved numerically. The mean-field ba- 
sis is found from the self-consistent HF equations: 

h(p)<p u = e v (j> v , (8) 

where 

h(p)=t + V(p) (9) 

is the single-particle HF Hamiltonian, t v are the single- 
particle energies, and p is the density matrix self- 
consistently determined by 

Pij = m^m- (io) 

The mean field potential is given by its matrix elements, 

v ij (p) = j2 v ikupik- (ii) 

kl 

In this conventional mean field formulation, the potential 
(fTTj) contains the direct and exchange contributions. The 
pairing effects, with strict particle number conservation, 
are contained in the superposition of the Slater deter- 
minants (2) used instead of the single HF determinant. 
The whole construction can be further improved by us- 
ing the exact but more complicated variational approach 
relating the single-particle basis to the full set of the co- 
efficients of the superposition (2). Such a possibility will 
be considered in the future. 

The HFP scheme of solution is the following: 

• Start with the spherical single-particle basis \ivm) 

• Choose in this basis the initial diagonal density ma- 
trix p corresponding to occupation numbers specific 
for prolate or oblate shapes (pairs with small or 
large \m\, respectively) 

• Solve the HF variational equation ((5|) and get the 
single-particle spectrum (<j) v ,e v ), in general corre- 
sponding to a deformed field 
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• Construct the "paired" class of many-body basis 
wave functions according to eq. © and calculate 
the matrix elements of the Hamiltonian H 

• Solve the variational equation ([7]) and obtain the 
ground state wave function 

• Calculate the next-step density matrix (fTQ|) 

• Repeat the procedure starting from the step three 
until convergence 

The converged results will certainly be a local mini- 
mum of eq. (2). Exploration of different starting choices 
in step 2 is needed to find a global minimum. In our 
study here we start with a spherical single-particle ba- 
sis (because the USDB interaction is so defined) but in 
principle any convenient axial basis could be used. In the 
end, the ground state energy can be found as the Hamil- 
tonian expectation value over the resulting ground state 
wave function |\&), 

£hfp = (*|#|*>. (12) 



III. RESULTS 

We performed calculations of ground state energies for 
all nuclei within the sd-shell region, from 17 to 39 Ca. 
Our results are summarized in Figs. 1-5. In Figs. 1 
and 2 we show the energy gain from HFP compared to 
HF, 

E CO rr = -EHF ~ -E?HFP . (13) 

Typical values are one to a few MeV. One observes the 
well-known odd-even staggering that is characteristic of 
pairing. In the conventional HFB approach the pairing 
correlation is zero for many of these nuclei, including 
cases such as 24 O where the spherical shell gap is too 
large, and cases such as 20 Ne, 24 Mg and 28 Si where the 
deformed shell gap is too large to support BCS type pair- 
ing. The HFP method gives some correlation energy for 
all sri-shell nuclei for which there are at least two ac- 
tive particles ( 9 < N < 19 or 9 < Z < 19 ). In the 
practical solution of the equations we find that many sd- 
shell nuclei have two or three energy minima. To have 
some confidence that we have found the lowest energy 
solution we start with several initial values of the den- 
sity matrix including those that are prolate and oblate 
deformed, spherical and random. 

In Figs. 3 and 4 we show the difference between the 
HFP result and the exact shell-model energy both ob- 
tained with the USDB Hamiltonian. For comparison of 
the methods, the full solution for the ground state of 28 Si 
must take into account 93,710 M-scheme Slater determi- 
nants. When projected onto good J — there are 9,216 
states, and when projected onto good J = 0, T = there 
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FIG. 1: Pairing correlation energy eq. (|13|) for all sri-shell 
nuclei. 
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FIG. 2: Four selected cuts on the correlation energies dis- 
played in Fig. 1 



are 839 states. The HFP method requires 92 determi- 
nants for protons and 92 for neutrons. The difference is 
clearly peaked at the N = Z line that can be explained 
by the proton-neutron pairing being not accounted for in 
the calculations. 

The HFP solution is very close to the exact solution 
around the edges of the sd shell (see Fig. 4). These nuclei 
are spherical and the HFP method is equivalent to the 
spherical exact- pairing method discussed in 0, Q ■ The 
largest deviation from exact is for nuclei near the middle 
of the sd shell. There are still pairing contributions for 
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FIG. 3: Energy difference between the exact results and our 
HFP model for all sd-shell nuclei. 
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FIG. 4: Energy difference between the exact results and our 
HFP model for four selected cuts through the sd shell. 



deformed nuclei, but the result is different from the naive 
expectation of just adding "spherical" contributions. For 
example, as shown in Fig. 2, the correlation energy is 
only about 400 keV for the deformed 20 Ne, compared to 
a total of about 3.4 MeV that would be obtained just 
from adding the 1.7 MeV correlation energies obtained 
for two neutrons and two protons in a spherical basis 
(e.g. 18 and 18 Ne). 

Finally in Fig. 5 we show the intrinsic quadrupole 
moment obtained for the lowest energy solutions for all 
sri-shell nuclei. One observes the well known region 



of strongly prolate nuclei near 24 Mg. 28 Si is the most 
strongly oblately deformed, and there is an island of weak 
oblate deformation around 31 Si. ft would be interesting 
to use our sri-shell sandbox to clarify the general question 
of why most nuclei are prolate deformed [8| , by exploring 
the HFP results with different (but realistic) Hamiltoni- 
ans. 
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FIG. 5: Intrinsic quadrupole moment Q, in units b 2 (b is the 
oscillator parameter), for the lowest energy configuration for 
all sd- shell nuclei. 



IV. CONCLUSION 



Obviously, the HFP is still far from adequate away 
from semi-magic nuclei. The angular momentum non- 
conservation is certainly a significant deficiency of the 
wave function, that when repaired will introduce addi- 
tional correlation energy. There are a number of ways 
that rotational correlation energies can be calculated, 
and we are optimistic that HFP wave functions can be 
used as a better starting point. 

Besides rotational energies, proton-neutron pairing ef- 
fects are omitted in our wave functions. As known 0], 
such pairing correlations are quite strong close to the 
N = Z line; they should be included at the next stage 
of development. Effective Hamiltonians for the HFB so- 
lution are explored in [lol ]. Finally, some improvement 
may follow from including the non-axial configurations 
with the pairing between more general time-conjugate 
orbitals (most probably, the mean field in 2 Mg is triax- 
ial MM)- 
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